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Abstract 

A few final comments on arXiv:1210.7548 are given to confute incorrect ar- 
guments claimed there. 

After our comment [T] on the arguments of ref. [2] , some of the same authors 
persist in some misleading arguments [3]. Here, as our last remarks on the topic, 
we comment a few points with particular emphasis to the arguments related to the 
extraction of the Dark Matter signal through the annual modulation signature and 
the role of the modeling of the background components in the single-hit counting 
rate. 

1 Extraction of the Dark Matter signal through the 
annual modulation signature 

Firstly, let us briefly recall a few arguments in order to focus the point. Two dif- 
ferent approaches are exploited in the field of Dark Matter (DM) direct detection. 
The first one aims at extracting the constant part (So) of the signal from the mea- 
sured counting rate, assuming a particular class of DM candidates. In experiments 
exploiting such an approach many kinds of uncertain subtractions /selections of the 
measured events are applied to select recoil- like candidates. 

The second approach, followed by the DAMA experiments, exploits a model- 
independent signature with very peculiar features: the DM annual modulation 
signature (see for example [HE])- in this case the experimental observable is not 
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So, but the modulation amplitude, S m , as a function of energy. This approach 
has several advantages; in particular, in this approach the only background of 
interest is that able to mimic the signature, i.e. able to account for the whole 
observed modulation amplitude and to simultaneously satisfy all the numerous 
specific peculiarities of the signature. No background of this sort has been found 
or suggested by anyone. 

Thus, the DM annual modulation model-independent approach does not re- 
quire any identification of So from the total single-hit counting rate, in order to 
establish the presence of DM particles in the galactic halo. 

So can be worked out only in a model dependent way, once S m has been 
determined. This procedure is performed by our collaboration by employing, 
within each specific framework, a maximum likelihood analysis which also takes 
into account the energy behaviour of each detector (see literature). 

In conclusion, the DM annual modulation signature allows one to over- 
come the large uncertainties associated to the exploitation of many data 
selections/subractions/statistical-discrimination procedures, to the modeling of 
surviving background in keV region and to the a priori assumption on the nature 
and interaction type of the DM particle(s). In particular, as already mentioned 
e.g. in [6], a precise modeling of background in the keV region counting rate is al- 
ways unlikely because e.g. of (i) the limitation of Monte-Carlo simulations at very 
low energies; (ii) the fact that often just upper limits for residual contaminants 
are available (and thus the real amount is unknown); (iii) the unknown location 
of each residual contaminant in each component of the set-up; (iv) the possible 
presence of non-standard contaminants, generally unaccounted; (v) etc.. 

A visual indication has, however, been given in ref. [7], where the cumulative 
energy spectrum over the 25 detectors has been reported up to 10 keV just as an 
example, showing that there was room for a sizeable constant part of the signal: 
namely So < 0.25 cpd/kg/keV in the [2,4] keV energy interval; this has been 
discussed by our collaboration many times in presentations at conferences and 
workshops. 

As a matter of fact, any attempt to extrapolate a modeling of the background 
in the keV single-hit counting rate can only demonstrate that room for Sq can 
exist, while it is not able to safely exclude a Dark Matter contribution. 

In conclusion, this explains why the conclusions derived from some fitting 
procedures used in ref. [3] are untenable (see also later). 

2 Other additional comments 

Other comments deserve to be reported in the following. 

Item regarding the so-called critique 1. We stressed in pQ and references 
therein that the electron capture of 40 K to the ground state of Ar - although its 
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branching ratio is not well known from the theoretical and experimental points of 
view - provides a small (10%) contribution to the total 40 K contribution at low 
energy in the single-hit counting rate. Moreover, the calculation of 40 K is not fully 
correctly performed in [21 [3], since e.g. the percentage of K shell electron capture is 
not taken into account. In particular, the probabilities of K shell electron capture 
to the first excited level (at 1461 keV) of 40 Ar (76.3% [8]) and to the ground state 
of Ar (87.9% [8]) must be included in the calculation. Furthermore that small 
(10%) contribution to the total K contribution at low energy in the single-hit 
counting rate has always been included by our collaboration in all the evaluations 

Item regarding the so-called critique 2. The authors in [3] seem to forget 
that the DAMA/LIBRA set-up is made of 25 detectors, each one with its own 
characteristics and energy spectrum. In the DAMA literature it was stated that 
"The analysis has given for the nai K content in the crystals values not exceeding 
about 20 ppb" [6]; as is evident, this was not an upper limit on nat K (with a 
relevant confidence level), but it is the maximum value among those measured in 
the 25 detectors. The average value (13 ppb) has been published in ref. [7] and 
discussed in many conferences. Therefore, there are in [21 [3] many statements 
("no data is presented to support this number, nor does the collaboration provide 
the uncertainty associated with it"; "the lack of details so far provided by the 
collaboration is unsettling.", ...) not justified and unfounded. 

Item regarding the so-called critique 3. A part from the arguments reported 
in Sect. 1 about the role of the modeling of the background components in the 
single-hit counting rate, we add here other specific comments. 

The authors of ref. [3] use different fitting functions (two segments + a gaus- 
sian) and different parameters with respect to our fit and their choice appears 
arbitrary. In particular, all the considerations in ref. [3] are also based on two 
questionable points: 

1. the assumption that the background in the low energy region is flat; 

2. the a priori belief that the level of this flat background in the [2,7] keV region 
can be arbitrarily fixed without taking into account any So contribution from 
the Dark Matter signaQ- 

As regards the first point, this assumption is justified in ref. [3] with "the 
universal feature of f3~ decays for small electron velocities", but this motivation 
is not sufficient. In fact, other components can contribute to the counting rate in 
the low energy region, and this is also evident in the given experimental spectrum, 

1 This assumption is always methodologically incorrect (and, in the particular case, also ex- 
perimentally contradicted by the measured observable S m and by the absence of processes able 
to mimic it). 
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that is far from being flat in the [2,10] keV region, as well as in the energy spectra 
published by other activities as shown e.g. in [7]. Moreover, it is odd that the 
authors (of ref. [3]) themselves violate their assumption by using a rising segment 
to describe the behaviour in the [7,10] keV low energy region. Different reasonable 
fits can provide slightly different upper limits for Sq, but the claim of a 0.85 
cpd/kg/keV flat background in the [2,7] keV energy region is completely arbitrary 
because it is not based on the knowledge of the background contributions but it is 
the result of a fitting procedure - among others - based on incorrect hypotheses. 

The authors of [3] conclude that a large modulation fraction (S m /So) is re- 
quired, but this is just the obvious consequence of their incorrect approach. 

Finally, just for completeness, let us note that in any case scenarios and Dark 
Matter candidates exist which can provide relatively large modulation fraction 
(see f.i. [TO]). 
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